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ABSTRACT 



We present and analyze the correlations between mid- infrared (MIR), far- 
infrared (FIR), total- infrared (TIR), Ha, and FUV luminosities for star-forming 
galaxies, composite galaxies and AGNs, based on a large sample of galaxies se- 
lected from the Spitzer SWIRE fields. The MIR luminosities of star-forming 
galaxies are well correlated with their Ha, TIR and FUV luminosities, and we 
re-scaled the MIR-derived SFR formulae according to the a bove correlations with 
differences less than 15%. We confirm the recent result by lCalzetti et al.l (120071 ) 
that the combined observed Ha + 24/im luminosities L(Ha Q b s + 24/rni) possess 
very tight correlation with the extinction-corrected Ha luminosities L(Ha_corr) 
for star-forming and even for dwarf galaxies, and show that the combined 
L(Ha Q bs+ 8/im[dust]) are also tightly correlated with L(Ha_corr) for the above 
sample galaxies. Among all the L(MIR)-L(FIR) correlations for star-forming 
galaxies, the L(24//m) vs. L(70yum) and L(8/xm[dust]) vs. L(160/xm) are the 
tightest and also nearly linear. The former could be related to young massive 
star formation, while the latter might be relevant to diffuse dust emissions heated 
by old stellar populations. Composite galaxies and AGNs have higher MIR-to- 
Ha/MIR-to-FUV luminosity ratios than star- forming galaxies, nevertheless their 
correlations among MIR, FIR and TIR luminosities are completely following 
those of star-forming galaxies. 
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stars: formation 
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Introduction 



Since the successful launch of Infrared Astronomical Satellite (IRAS) in 1983, more 
advan ced space-based infr ared telescopes had been launche d, such as In f rared Space Observatory 
(ISO; iKessler et al.lll996l ) and Spitzer Space Telescope (IWerner et al.l 120041 ). to help as to 
study the MIR to FIR emission properties of galaxies. MIR to FIR emissions are crucial for 
quantifying extragalactic star formation activities, since infrared emissions are arising from 
the dust re-radiation of UV photons emitted by young massive stars. FIR luminosity has 
been proved to be a good star f o rmation rate (SFR) indicat or of galaxies by IRAS obser- 



vations (e.g.. iHunter et al.lll986l ; iLehnert fc Heckman II1996T ). However, the low sensitivity 



and spatial resolution of FIR observations obstructed us from applying it to most galaxies 
at intermediate and high redshifts, even though there are some degree of improvements in 
the ISO and Spitzer era. Optical/UV SFR tracers, such as Ha-line emission and FUV con- 
tinuum luminosities, are important because they are still the only way to measure SFR in 
regions of low e xtinction such as fa int, dwarf galaxies and especially extended UV disk (XUV- 
disk) galaxies (IThilker et al.l 120071 ). But these tracers usually suffer from serious extinction 
effects for norm al and especially infrared-bright galaxies which are difficult to be corrected 
(|Jonsson!l2004l ). MIR emissions of normal galaxies are dominated by du st continuum from 
very small grains ( VSGs), together with some broad emission features (IGillett et al.lll973l ; 



Willner et al.ll 19771) which were realized to b e from polycyclic aromatic hydrocarbons (PAHs; 
Leger fc Puget Ill984l ; iPuget fc Legerlll989l ). Early analysis of correlation s between MIR hi 



minosity (der i ved f r om ISO observations) and o t her SFR tracers ( e.g., lElbaz et al.l 12002 
Roussel et al.l l200ll ; iForster Schreiber et al.l 12004 ; iFlores et al.l |2004j) have shown that the 
MIR luminosity could be used to estimate SFRs of galaxies. 

Mid-Infrared observations by Spitzer with higher sensitivity and better angular reso- 
lution than ISO, provide a new opportunity to study young stellar populations and star 
formation proces ses in galaxies. The 8/im band of Spitzer Infrared Array Camera (IRAC; 
Fazio et al.l 12004 ) is designed to cover the 7.7/im P AH feature: whil st the 24/im band of 
Multiband Imaging Photometer for Spitzer (MIPS; iRieke et al.l 12004 ) just covers the dust 
continuum from VSGs for local galaxies, avoiding the contamina t ion o f most of the MIR 
emission or absorption lines. The pioneering work by IWu et al.l ( 120051 ) (hereafter Wu05) 
has shown that, for star-forming galaxies, both Spitzer IRAC 8/zm(dust) and MIPS 24/im 
luminosities are strongly correlated with the 1.4 GHz radio and Ha luminosities. Then the 
new formulae for estimating SFRs of galaxies were derived fr om MIR luminositie s based on 
above correlations. Detailed studies of nearby galaxies M 51 (jCalzetti et al.ll2005l ) and M 81 
( iPerez-Gonzalez et al.ll2006l ) revealed that the 24/im luminosities of star- forming regions in 
the above two galaxies were well correlated with the extinction-corrected Paa luminosities. 
Based on the 8/im and 24/im data for HII regions in 33 nearby SINGS (Spitzer Infrared 
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Nearby Galaxies Survey; iKennicutt et al.l 120031 ) galaxies, ICalzetti et al.l (120071 ) found that 
both 8/xm 
They and 



dust) and 24iim lumin osities are correlated with L(Paa_corr), but non-linearly. 



Kennicutt et al.l (120071 ) also demonstrated that the combination of observed Ha 



(uncorrected for dust extinctions) and 24/im luminosities could be one of the best SFR 
indicator for HII regions in these galaxies. 

Besides being heated by young massive stars, the energy source of MIR and FIR emis- 
sions could also be the active galactic nuclei (AGN) in the galactic center, since AGNs are 
also powerful U V emitters. Weak AGN-hosting galaxies are very popular in the local universe 
(IHo et al.lll997l ). thus the possible AGN contribution to MIR and FIR emissions could have 
significant impa ct on the estima t ion of cosmic SFR density and the understanding of its evo- 



lution (see, e.g, iBell et al.ll2005l ; iPerez-Gonzalez et al.ll2006l ; lElbaz et al.l 120071 ; iDaddi et al. 



20071 ). Previous works showed that powe rful AGNs can streng t hen MIR VSG emissions 



but weaken or eve n avoid PAH emissions (jWeedman et al.ll2005l ; ISiebenmorgen et al.l 12004 



Verma et al.ll2005l ). Weak AGNs are f o und to have M IR properties different from that of 



star- forming galaxies (jWen et al.l 12007 ; 



could be related to the AGN activity (jWu et al 



Li et al. 



2007J), and the degree of such differences 
20071 ). 



The Spitzer Wide-area Infrared Extragalactic Survey (SWIRE; lLonsdale et al.ll2003l ) is 
the largest extragalactic survey program among the six Spitzer cycle- 1 Legacy Programs. 
With a total field of ~49 de g 2 , it is m u ch lar ger than the Spitzer First Look Survey (FLS) 
field (~3.7 deg 2 ) studied by IWu et al.l (120051 ). In this paper, we perform statistical analysis 
on the correlations between MIR, FIR, Ha, and FUV luminosities for star-forming galaxies, 
composite galaxies (which were affected by both central nuclear and star formation activi- 
ties), and AGNs selected from three northern SWIRE fields. 

The structure of the paper is as follows. We describe the construction of our sample 
and the estimation of multi-wavelength luminosities in §2 and §3. The major results on 
correlation analysis are presented in §4. Some discussions and a summary of this work are 
given in §5 and §6. Throughout this paper, we adopted a ACDM cosmology with f2 m = 0.3, 
Q A = 0.7 and H = 70km S - 1 Mp C - 1 . 



2. Sample and Data Reduction 

2.1. SDSS Sample 

Th e optical spectra l sample galaxies are selected from the S loan Digital Sky Survey 
(SDSS; lYork et al.l 120001 ) main galaxy sample ( Strauss et al.l 120021 ) . with r-band Petrosian 
magnitudes less than 17.77 mag. The SDSS main galaxy sample covers three northern 



-4- 



Spitzer SWIRE fields: most of the Lockman Hole (center: 14 41 00, +59 25 00) and ELAIS- 
N2 (center: 16 36 48, +41 01 45), and about 1/3 of ELAIS-N1 (center: 16 11 00, +55 00 
00). The overlay regions used in this paper is about 15 deg 2 . A total of 821 galaxies were 
selected with available Ha and H/5 emission line fluxes (> 5a), based on the emission line 
catalog of SDSS DR40. 



2.2. MIR Sample 



The IRAC four bands (3.6, 4.5, 5.8 and 8.0 /xm) images were mosaicked from the Basic 
Calibrated Data (BCD; acquired from the Spitzer S ciences Center) a f ter flat-field correc - 



tions, dark subtraction, linearity a nd flux calibrations (IFazio et al. 



with the final pixel scale of 0.6" (jWu et aUl2005t ICao & Wu 



2004 



2007 



Huang et al l2004h 



Cao et al|l2008j ). The 



MIPS 24/xm images were mosaicked as the same way but with the pixel scale of 1.225 



( IWen et al.ll2007l : ICao & Wu 120071: IWu et al 



20071 ). Based on the catalogs of the Two Mi- 



cron All Sky Survey (2MASS; ICutri et al.ll2003l ). we obtained the accuracy of the astrometric 
calibration around 0.1" in all the five bands. 

We obtained the auto and aperture magnitudes for all the sources in the three SWIRE 
fields detected by SExtractor (IBertin &: Arnoutsl Il996l ). The auto magnitudes were used 
to denote the total luminosities of detected sources. Hereafter, the auto magnitudes were 
represented by the total magnitudes for these five MIR bands. Aperture photometries with 
3" and 6" diameters were adopted for the IRAC four bands and MIPS 24/im band images, 
respectively. 

The aperture corrections for the IRAC four band photometries (corrected to an aper- 
ture of 24") are: —0.66, —0.64, —0.80 and —1.00, respectively. We re-calibrate d the IRAC 
apertu re magnitudes with the aid of 2MASS photometric catalog, according to lLacy et al. 
f l2005h and lEisenhardt et all 120041 ). Firstly, we selected blue 2 MASS field stars (F7 dwarfs) 
with J — K color less than 0.3. To ensure accurate photometry and avoid saturated sources 
in IRAC 3.6/zm images, the K-band magnitudes were restricted within 10 to 14 mags; and 
then, the measured AB magnitudes in each IRAC band were converted to Vega magnitudes 
by adding factors of 2.79, 3.26, 3.73 and 4.40 ft res pectively; afterwards, we compared the 
four K — IRAC colors with 1 Expected Value'' from lBessel fc Brettl (119881 ) and adopted the 
averaged values of 0.04, -0.04, -0.03 and 0.07 (based on Table [T]) as the corrections for each 



1 http://www. mpa-garching.mpg.de/SDSS/DR4/; see also 
2 http:/ /spider. ipac.caltech.edu/staff/gillian/cal. html 



Tremonti et al. 



(120041 ). 
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band. The final flux calibration uncertainties of all the four IRAC bands are less than 0.08 
mag. The 6" aperture photometry of MIPS 24/mi band was corrected to an aperture of 30" 
, with an aperture correction value of —1.21. Following the procedure used by the MIPS 
instrument team to derive calibration factors from standard star observations, we corrected 
an additional calibration factor of 1.15. A factor of 0.96 was divided in order for vF v equat- 
i ng constant scale d The final MIPS 24 /im flux has a calibration uncertainty less than 10% 
(kieke et al.ll2004f ). 



To construct the MIR galaxy sample, we cross-matched these MIR sources with the 
above SDSS sample galaxies with a radius of 2". Only 15 out of the 821 SDSS galaxy are 
mismatched. Most of these sources are edge-on galaxies or have bright neighbors. The 
completeness of the MIR sample is ~98%, since the sensitivity of Spitzer MIR bands is 
much deeper than that of SDSS spectroscopic observations. For the purpose of this work, we 
need both total and aperture magnitudes in at least four MIR bands: IRAC 3.6, 5.8, 8 /im, 
and MIPS 24/im. Therefore, another 201 galaxies were rejected because of their locations at 
the edge of overlay regions of SDSS and/or Spitzer images. The final number of our MIR 
sample galaxies is 605. 



2.3. FIR Sample 

The Spitzer MIPS 70 and 160/im band photometry catalogs of SWIRE northern fields 
were derived from the SWIRE data release 3 @. The 70/im catalog was cross-matched with 
the 160/zm catalog with a radius of 11", and then cross-matched with the MIR sample with 
a radius of 4". Since we need to calculate the total infrared (TIR) luminosity, only the 
sources with detectable fluxes in all three MIPS bands were selected. The final FIR sample 
includes 197 galaxies. The point-response function (PRF) magnitudes were adopted as total 
magnitudes in this work, and the absolute flux uncertainties are 20% and 30% for the two 
bands, respectively B 



3 http: / /swire. ipac.caltech.edu/swire/astronomers/publications/SWIRE2_doc_083105.pdf 
4 http:/ /swire. ipac.caltech.edu/swire/astronomers/data_access. html 

5 http:/ /swire. ipac.caltech.edu/swire/astronomers/publications/SWIRE2_doc_083105.pdf 
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2.4. UV Sample 



The Galaxy Evolution Explorer (G 'ALEX; iMartin et al.ll2005l ) General Data Release 
2 (GR2) and General Data Release 3 (GR3) have been released to public in spring 2006 
and summer 2007. Two UV bands were used: the far-UV (FUV) band covering ~ 135-175 
nm and the near-UV (NUV) band covering ~ 175-275 nm. The UV catalogs were taken 
from pipeline-processed GALEX Deep Imaging Survey (DIS) (hiab~25; IMartin et al.ll2005l ) 
observations of three northern SWIRE fields. The sources with calibrated Kron magnitudes 
in both NUV and FUV bands were selected. After cross-matched wit h the MIR sample w ith 
radius of 4", 421 galaxies were left. We need to use the formula from lTreyer et al.l (120071 ) to 
compute the intrinsic extinction of the sources, but this formula is not suitable for the red 
sequence galaxies, which were defined with the NUV — r color < 4 (at redshift 0.1). Finally, 
359 galaxies were included in our UV sample. The calibration accuracy of the measurements 
is about 3% H 



2.5. Spectral Classification 



The optical spectral classifications of galaxies were made by adopting the traditional 



BPT diagnostic diagram: [Nil] /Ha versus [OIII]/H/3 (Baldwin et alll981 



Veilleux fc Osterbrock 



19871 ) , as shown in Figure [H The dashed curve is from iKauffmann et al.l (120031 ) (hereafter 
Ka03 line) and the dotted curve is from iKewlev et all (bpOlh (hereafter KeOl line). The 
objects located below the Ka03 line were classified as star-forming galaxies; those above 
the KeOl line were classified as narrow-line AGNs (note we did not distinguish Seyfert 2s 
from LINERs (iHeckmanl Il980l ) here); thos e between the aboy e two lines were classified 
as composite (starburst + AGN) galaxies ( Kewlev et al.l 120061) . Composite galaxies are 
similar to mixture-type galaxies defined by 



WuetaL 



(jl998j ) 



which were based on three 
diagnostic diagrams (adding [SII]/Ha vs. [OIII]/H/3 and [OI] /Ha vs. [OIII]/H/3) from 
Veilleux fc Osterbrock I (119871 ): if the classification of a source in the three diagnostic dia- 
grams is different (e.g., it was shown in one diagram as a star-forming galaxy, but could 
be classified as a AGN by the others), then it was categorized as a mixture- type galaxy 
( IWu et al.lll998l ). The numbers on spectral classifications of our MIR, FIR and UV samples 
are shown in Table [2J 



6 http: / /galex.stsci.edu/GR2 / 

7 lit t p : / / galexgi . gsfc . nasa. gov /docs/ galex /Cycle3/ 
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2.6. Sample Distribution 



The distributions of the SDSS B-band absolute magnitude (Mb), redshift, u — r color, 
and equivalent width of Ha line emission (EW[Ha]) for our MIR, FIR and UV sample 
galaxies are shown in Figure [2j Most of the galaxies have Mb brighter than —18 mag and 
redshift less than 0.2. Only a few of them hav e magnitudes lower than —18 mag, which 
belong to dwarf galaxies (IThuan fc Martin 1981). Here, Mr were calculated from the SDSS 
g and r-band magnitudes according to lSmith et al.l (120021 ). The three galaxy samples do not 
show obvious differences in all of the above distributions. 

Figure [3] shows the above four distributions for three different spectral type galaxies 
(star- forming galaxies, composite galaxies, and AGNs). No obvious difference was found 
in the distributions of either Mb or redshift. However, the distributions of u — r color and 
EW(Hct) for composite galaxies and AGNs are strikingly different from those for star-forming 
galaxies. This result cou ld be due to the fact that AGNs tend to reside in earlier type galaxies 
( IKauffmann et ajj [2003| ; see §5.4 for detailed discussions). 



Estimation of Luminosities 



3.1. MIR luminosities 



To obtain MIR dust emission luminosities, the MIR continu um contributed by star light 
should be firstly estimated and then removed. A factor of 0.232 (IHelou et al.ll2004j ) was used 
to scale the stellar continuum of 3.6/xm to that of 8/xm, with the assumption that the entire 
IRAC 3.6/xm band emission is from stellar emission, and based on the Starburst99 synthesis 
model (ILeitherer et al.lll999l ). assuming solar metallicity and a Salpeter initial mass function 
(IMF) between 0.1 and 120 M Q . The scaling factor used here is not quite different from 
that of 0.26 by Wu05, which was obtained based on statistics of elliptical galaxies in the 
FLS field. Following Wu05, we denoted 8/zm(dust) as the 8/xm dust emission. The star light 
contribution to the MIPS 24/im band can be neglected since it decreases quickly with the 
increase of wavelength. 

Since IRAC 8/im band cover the complex PAH emission range, a more accurate spectral 
energy distribution(SED) mo del is needed t o obt ain a reliable MIR K-correction. Here we 
adopted a series of SEDs from lHuang et al.l (|2007|)'s two-component model which is a linear 
combination of an old 'early-type' stellar population and a 'late-type' spiral disk population. 
The MIR color [5.8]- [8.0] was used to select the best SED model, and then the K-correction 
value was obtained. A template SED of a normal HII galaxy NGC 3351 (from SINGS; 
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Kennicutt et al.ll2003l ) was used to perform K-correction for the 24/zm band for all the sample 
galaxies. 



3.2. FIR and TIR luminosities 

The K-correction of FIR 7 and 160/xm emissions was b ased on the best fitted SED 
model of starburst galaxy M82 (jSiebenmorgen fc Kriigell 120071 ). The correction factors are 
less than 12% and 30% for the ab ove two bands. The TIR luminosity (3-1100yum) was 



calculated using the Equation (4) of iDale fc Helou I (120021 ): 

L[TIR] = 1.559 x ^[24/1 m] + 0.7686 x uL u [70fxm\ + 1.347 x z/L„[160//m] 



3.3. Ha luminosity 

The optical spectra of SDSS were taken with 3" diameter fibers, thus the measured 
Ha emission fluxes were only from the central regions for most low redshif t galaxies and 



the co rresponding aperture correction is needed. We adopted the method by iHopkins et al. 



(120031 ) and used their Equation (A2) to obtain the Ha luminosity of the whole galaxy: 

L[Ha) = AnDl S Ra lO" - 4 ^--^-) } (2) 

Here, D L is the luminosity distance; Sua is flux of Ha emission; while r Petro and rgb er are the 
r-band Petrosian and fiber magnitudes, respectively. To match the 6" aperture photometry 
in 24/iin band, we also calculated the Ha luminosities in aperture of 6" in a similar way, 
i.e., substituting the r-band 6" aperture magnitudes for rp etro in the above equation. The 
r-band 6" aperture magnitude was derived from the r-band radial surface brightness profile 
of the SDSS catalog Jstoughton et al.ll2002h . 



The optical Ha emissions suffer from the dust extinctions from both Milky Way and 
the host galaxy. Th e foreground Galactic extinction was first corrected by assuming the 
Cardelli et al.l (119891 ) 's extinction curve and Ry=3.1. Then the intrinsic extinction correction 
was performed, base d on the colo r excess E(B-V) which was acquired from the Balmer 
decrement F Hq /F H /3 (ICalzettill200ll ). 



3.4. FUV luminosity 



The K-correction of UV band was based on the IDL code supplied by Bl anton (version 
tag v4_l_4). The method and the SEDs used by this code were described by iBlanton et al. 
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( 120031 ) and lBlanton &: Roweia (120071 ). Since dust will heavily obscure the UV emission, it is 
important to correct f or the UV exti n ction . The Galactic extinction was determined by the 
par ametrized curve of ICardelli et al.l (119891 ). with the adopt ed conversion fac tors of 7.9 and 
8.0 JGil de Paz et al.ll2007h for FUV and NUV, respectively fauat et alJ l 2005h. The intrinsic 



UV extinctions a re usually computed by the TIR-to-UV ratio method ( Kong et al 



Buat et al. 



2005. 



20051 ; ISeibert et al.l 120051 ; ICortese et al.l 120061 ; IZamojski et al. 



2007 



2004 



Salim et al. 



20071 ). However, the number of our UV sample galaxies with FIR data is too small 



to make statistical analysis. Alternatively, using a large sample of ga l axies from SDSS 
spectroscopic catalog and GALEX Medium Imaging Survey, iTreyer et al.l (120071 ) derived an 
empirical correlation to estimate the intrinsic UV extinction based on the FUV — NUV 
color at redshift ~ 0.1 (their equation (9)): 



A-fuv = 4.05 x (m FUV - m NUV ) - 0.18. 



(3) 



We use this equation to calculate the intrinsic UV extinctions for our UV sample galaxies. 



4. Correlation Analysis 
4.1. L(MIR) vs. L(Har) 

Total luminosity 

Firstly, we present the correlation between total MIR luminosities and total Ha lumi- 
nosities for star-forming galaxies. The total Ha luminosities are the aperture-corrected values 
(see §3.3). Total 24/xm and 8/im(dust) luminosities are plotted against total Ha luminosi- 
ties in Figure H] (a) and (b). After excluding the dwarf galaxies (open circles), star- forming 
galaxies show good correlations between total L(MIR) and L(Ha). Using the two-variable 
regression, we obtain the best nonlinear fits (log 10 (y) =a+blog 10 (a;); solid lines in Figure @] 
(a) and (b)): 

Log 10 vL v [2Aum\ total = (-0.06 ± 0.07) + (1.18 ± 0.02) x Log 10 L[H a ] total , (4) 

Log w vL„[8.0um(dust)] total = (1.14 ± 0.07) + (1.07 ± 0.03) x Log w L[H a ] total . (5) 

Comparing with the nonlinear fits for star-forming galaxies in the Spitzer FLS field obtained 
by Wu05 (dashed lines in Figure 0] (a) and (b)), we find that the best nonlinear fit of total 
24/im vs. Ha luminosities in this work is well consistent with the corresponding result by 
Wu05. However, the slope of our nonlinear fit of total 8/mi(dust) vs. Ha is a little shallower 
than that of Wu05, and this discrepancy could be due to the 8/im-weak HII galaxies (located 
at bottom-left of Figure H] (b)) which tend to have lower 8/mi(dust)-to-Ha luminosity ratios. 
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Similar to Wu05, we can derive MIR SFR formulae based on the best linear fits ( 
\og 10 (y) =c+ log 10 (a;); dotted li nes in Figu r e H] (a ) and (b), see also Table [3] for the fitting 



parameters) and Formula (2) of iKennicuttl (119981 ) 



SFR(M & yr~ l ) = 7.9 x W' 42 LlH^ergs- 1 ), (6) 

Here, solar abundances and the Salpeter IMF (0.1-100 M ) were used in deriving the cali- 
bration factor. So the derived formulae according to the above correlations are: 

SFR^ m {M Q yr ) = ^ - ^ , (7) 



SFR 8 , m(dmt) (M Q yr ) = L5g x 1Q9 , 



J tota/ / \ 
. (8) 

The fitting coefficients are consistent with those derived by Wu05 with differences of ~11%. 



Aperture luminosity 

In the above analysis, the aperture effect was corrected and then the total Ha lumi- 
nosities were obtained based on the assumption that the outer parts of galaxies have the 
same Ha-to-r luminosity ratios to the central regions. In fact, this could not be true for 
many reasons, e.g., for central starburst galaxies, it would overestimate the total Ha fluxes. 
To evaluate such an effect, we re-analyzed above correlations only for the central regions 
of star- forming galaxies. Figure H] (c) and (d) show the correlations between L(MIR) and 
L(Ha) in aperture of 6" at 24 fim and of 3" at 8/zm. The best nonlinear fits are: 

Log l0 uL u [2A^m] aper = (0.63 ± 0.05) + (1.10 ± 0.02) x Log l0 L[H a } aper , (9) 

Log 10 uL u [8.0fim{dust)} aper = (1.74 ± 0.05) + (1.02 ± 0.02) x Log 10 L[H a } aper . (10) 

Obviously, the scatters shown in Panels (c) and (d) (the standard deviation of the fitting 
residuals are 0.15 and 0.14 respectively) are smaller than those in Panels (a) and (b) (the 
standard deviation of the fitting residuals are 0.16 and 0.18 respectively) of Figure HI es- 
pecially for the correlation between 8//m(dust) vs. Ha luminosities (see also the fitting 
coefficients summarized in Table [3]). Compared with the best nonlinear fits for total L(MIR) 
vs. L(Ha) (dashed lines in FigureH](c) and (d)), the best nonlinear fits for aperture L(MIR) 
vs. L(Ha) (solid lines in Figure H](c) and (d)) are about 0.1 dex shifted upward. This result 
indicates that the MIR-to-Ha luminosity ratios of the central regions are relatively higher 
than those of the whole galaxies. Panels (c) and (d) ( "aperture" correlations) also show that 
the deviations of dwarf galaxies from normal star-forming galaxies are much less than those 
shown in Panels (a) and (b) ("total" correlations). 
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4.2. L(Ha obs + MIR) vs. L(Hot) 



Calzetti et al.l (120071 ) showed that there existed a tight correlation between the linear 
combination of the observed Ha (Ha b s ) and 24/zm luminosities and the extinction-corrected 
Paa luminosities for star-forming knots in SINGS galaxies. This result can be interpreted 
as: the 24/xm (and also 8/im[dust]) and observed Ha emissions trace the dust-obscured and 
un-obscured s tar formation, respecti vely. This correlation could be presented as equation 
illustrated by lKennicutt et al.l (120071 ) (their equation 4): 

L[H a ] = L [H a _ obs ] +axvL u [24/im] . (11) 

Here, the left side of this equation is the extinction-corrected Ha luminosities. 



Based on these star-forming galaxies, we can obtain the scaling constant a (jKennicutt et al. 



20071 ) values of 0.022 and 0.010 for the combined Ha Q b s + 24/im and Ha Q b s + 8/im(dust) 
luminosities. The a value of 0.022 we obtained here i s sligh tly smaller than that of 0.031 
for star-forming regions in galaxies by ICalzetti et al.l (120071 ). and we speculate that this 
difference is caused by the dust heated by the interstellar radiation field (mainly from old 
stars) which could enhance the MIR emissions for the entire galaxies. In addition, the tight 
correlations between the dust-obscured (absorbed) Ha luminosities and both of the 24/im 
and 8/im(dust) luminosities also verify the fact that the dominant sources of MIR emissions 
are from the dust-obscured star formation because the Ha emission directly traces the HII 
region surrounded young massive stellar population. Therefore, the linear combinations of 
the Ha bs and MIR luminosities would include almost all the star formation activities in 
galaxies, and they could be used as good SFR indicators for normal star-forming galaxies. 

With the derived constant a, the tight correlations between combined Ha obs + MIR 
and extinction-corrected Ha luminosities (see Table [3]) are shown in Figure [5j The dwarf 
galaxies also follow the correlations for normal star-forming galaxies. The best nonlinear fits 
are shown as: 

Log 10 (L[H a _ obs ] + 0.022 vL u [2^f£m]) = (0.29 ± 0.04) + (0.96 ± 0.01) x Log 10 L[H a ). (12) 
Log 10 (L[H a _ obs ] + 0.010 uL u [8fim]) = (0.86 ± 0.04) + (0.89 ± 0.01) x Log w L[H a }. (13) 

Thus, the combined L(Ha D b s + 24/im) almost linearly correlate with extinction-corrected 
Ha luminosities. 



4.3. L(MIR) vs. L(FIR) 



FIR luminosity was thought to be a good SFR tracer by previous IRAS and ISO 
observations. Figure [6] shows the correlations between MIR (24/aul and 8/mi[dust]) and FIR 
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(70/im and 160/im) luminosities. Due to the relatively lower sensitivities of Spitzer MIPS 
70 and 160/im observations, only 141 star-forming galaxies have both 70/zm and 160/im 
detections. The best nonlinear and linear fits are shown as the solid and dotted lines in 
Figure and the fitting parameters are listed in Table [31 Among the above four correlations, 
the L(24/irn)-L(70/im) (Figure [6](a)) and L(8/xm[dust])-L(160/xm) (Figure EJ^d)) are shown 
to be much tighter than the rest two correlations (Figure [6](b) and (c)). The former two 
correlations have scatters with standard deviations of ~0.07, while those of the latter two 
are about 0.10. Furthermore, from Table [31 it is clear that the slopes of best nonlinear fits of 
L(24/xm)-L(70/xm) and L(8/im[dust])-L(160/im) correlations are both ~1.04, which indicate 
that these two correlations are almost linear. 



4.4. L(MIR) & L(FIR) vs. L(TIR) 



The correlations between MIR & FIR with TIR (3-1100/xm) luminosities are shown in 
Figure [71 The fitting parameters and correlation coefficients are listed in Table [31 All of the 
four correlations are tight, with the standard deviations of 0.06 to 0.09 and the Spearman 
Rank-order correlation analysis coefficients of ~0.97. The L(160/xm)-L(TIR) correlation is 
the tightest among the four, and the 8/mi(dust) and 160/im luminosities are almost linearly 
correlated with the TIR luminosities. 



Based on these correlations 
nosities. 



we re-calibrated the MIR SFR formulae from TIR lumi- 



Kennicuttl (119981 ) presented the equation (Equation (4)) to calculate SFR based on 
IR luminosity: 

SFRiMayr- 1 ) = 4.5 x IQ- U L[F I Rings' 1 ), (14) 



The calibration factor was derived by applying the models of iLeitherer fc Heckman I (119951 ) 
for continuous bursts of age 10-100 Myr and adopting the same Salpeter IMF. Here, Lpm 
refers to the IR luminosity integrated over the full-, mid-, and far-IR spectrum (8-1000 /im). 
The TIR (3-1 100/im) luminosity c an also be derived from IRA S three bands (25,60,100/xm) 
luminosities (IDale fc Helou 1 120021 ). and iTakeuchi et al.l (120051 ) found that the L(TIR) es- 
timated in this way w ere very similar to the L (IR) (8-1000/iin; calculated by using the 
equations presented by ISanders fc Mirabel Ill996l ). Therefore, we adopted above Equation 
and used the TIR luminosities to calibrate SFRs. The derived formulae are: 

z/Lj,[24/xm] 



S FR 2 4fim ( M Q yr ' 



SFR, 



8/^m(dust) 



(M & yr- 



5.62 x 1O 8 L ' 
uL v [8fj l m(dust)] 



1.17 x 1O 9 L ' 

The fitting coefficients are different from those in Equations (6) and (7) by 



(15) 



(16) 



20-25%. 
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Kennicuttl ( 119981 ) pointed out that the calibration factor of 4.5 in their Equation (4) 
was not valid for quiescent, normal galaxies, since the lower optical depth of dust could 
increase this coefficient, while the dust heating from old stars could decrease it. The old 
stellar population could heat the interstel lar dust and be re-radiated at longer wavelength 
as cool infrared cirrus emissions (see, e.g.. lBeichmanlll987l ; IChoill2006l ). Therefore, the SFR 



derived from TIR lumi n ositie s would be overestimated if such contributions to L(TIR) can 
not be neglected. iBelll (120031 ) illustrated that the fraction of contribution to L(TIR) from 
diffuse components was ~32% for normal star-forming galaxies, while this proportion is much 
smaller for Luminous Infrared Galaxies (LIRGs) with L(TIR) > lO n L . This interpretation 
could explain the differences between fitting coefficients of SFRs calibrated by TIR and Ha 
luminosities. 

Considering that most of our sample galaxies have infrared luminosities less than 10 11 
Lfp, w e adopt an old stellar population fraction of ~30 % independent of luminosity (IHopkins et al 
20031 ). Then we obtain old stellar population corrected SFR formulae for star- forming galax- 
ies (with L(TIR) < lO n L ): 

uL u [24/xm] 



SFR 24tim (M Q yr- 



SFR 



8/^m(dust) 



(M Q yr- 



7.79 x 1O 8 L ' 
uL u [8fim(dust)} 



(17) 



'11 



1.70 x 1O 9 L 

The fitting coefficients are consistent with those in Equations (6) and (7), with much smaller 
differences of ~8%. 



4.5. L(MIR) vs. L(FUV) 

Figure [8] shows the 24/xm and 8/im(dust) vs. extinction-corrected FUV luminosities (see 
§3.4) for 248 star-forming galaxies. Although scatters are large (with standard deviations of 
~0.22), there still exist correlations between L(MIR) and L(FUV) with the Spearman Rank- 
order correlation analysis coefficients of ~0.86, 0.8 4, respectively. T he conversion between 



the UV luminosity and SFR has been presented by iKennicuttl (119981 ) (Equation (1)): 



SFR(M Q yr- 1 ) = 1.4 x 10~ 28 L u (erg S - l Hz~ l ), (19) 

The calibration factor was derived by the assumption of continuous star formation and 
adopting the same Salpeter IMF. Based on this formula, we derive SFRs from the FUV 
luminosities and the linear correlation coefficients listed in Table [3j 

W')^, (20) 
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SFR 



8^tm(dust) 



(M e yr 



-l^ 



uL u [8fim(dust)} 



(21) 



0.93 x 1O 9 L 

The fitting coefficients are nearly half of those in Equations (6) and (7), but note the FUV 
luminosities only cover a range of ~2 orders of magnitude. 



However, iKennicuttl (119981) 's formula is only su it for continuous star formation over time 
scales of 10 8 year or longer. iBurgarella et al.l (120051 ) presented a form ula (their Equation (2)) 
to estimate SFRs for non-starburst galaxies whose stellar birthrate b (IKennicutt et al.111 994) 
are les s than 1.0. Based on the relationship between b and EW(Ha) given by IKennicutt et al. 
(119941 ). we selected 187 star- forming galaxies whose EW(Ha) are less than 30 (correspond 
to b < 1) to re-calibrate SFRs: 



SFR 2ilim (M & yr- 1 ) 



uL u [24:fim\ 
7.04 x 1O 8 L ' 



(22) 
(23) 



Similar to the analysis in §4.2, Figure [9] presents the correlations between the extinction- 
corrected FUV luminosities and combined observed FUV + MIR luminosities (L[FUV Q b s + 
MIR]), with derived coefficients a of 6.31 and 3.02 for FUV + 24/mi and 8/im, respectively. 
But the scatters of these correlations (~1.65) are larger than that of L(Ha) vs. L(Ha; bs + 
MIR). 



4.6. Correlations for Composite galaxies and AGNs 

Besides the star-forming galaxies, there are also many composite galaxies and AGNs in 
our sample. The correlations between the total & aperture MIR and Ha luminosities for 
composite galaxies and AGNs are shown in Figure [101 The parameters of the best nonlinear 
and linear fits for composite galaxies are listed in Table HI The L(MIR) are correlated with 
L(Hct) for both composite galaxies and AGNs, but the scatters are larger (with standard 
deviations of 0.15-0.21) and the correlation slopes are shallower than those of star-forming 
galaxies. Furthermore, the MIR-to-Ha luminosity ratios of both composite galaxies and 
AGNs are higher than that of star- forming galaxies, especially for the L(24/im)-to-L(Ha) 
ratios. Similarly, both composite galaxies and AGNs are also shown to have relatively higher 
MIR-to-FUV luminosity ratios (see Figure [S]). 

Figures E] and [H show the correlations of L(MIR) vs. L(FIR) and L(MIR) & L(FIR) 
vs. L(TIR) for composite galaxies and AGNs. The best nonlinear and linear fits for com- 
posite galaxies are similar to those for star-forming galaxies but with larger scatters. Just 
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the same as normal star-forming galaxies, for composite galaxies the L(24//m)-L(70//m) 
and L(8yum[dust])-L(160/im) correlations are tighter than those of L(24/xm)-L(160/xm) and 
L(8/xm[dust])-L(70/im); the best nonlinear fits of the former two correlations are almost lin- 
ear; and both 8/xm(dust) and 160/im luminosities are almost linearly correlated with TIR 
luminosities. 



5. Discussion 



5.1. Discrepancies in L(MIR)-L(Ha) Correlations 



The aperture correction to Ha luminosities (see §3.3) were performed based on the 
assumption that the distributions of Ha emission in galaxies are similar to those of the 
continuum (SDSS-r) emission. All of our sample galaxies have both detectable Ha and B./3 
emission lines above 5a significance level, thus it could bias towards galaxies with stronger 
(circum-)nuclear star formation activities. In the outer parts of these types of galaxies, star 
formation activity would be much weaker and the optical continuum emissions could be 
dominated by old stellar populations. In such cases, we would overestimate the total Ha 
luminosity if performing corrections from SDSS r-band magnitude. The aperture correction 
effect discussed above could explain the shifts between best fits of total and aperture MIR- 
to-Ha luminosity correlations shown in Figures H] and described in §4.1. 

Besides the aperture correction effect discussed above, differential dust obscuration could 
be another reason for such discrepancies. The central region of gal axies is more dusty 
than their outer parts (e.g., iPopescu et al.l 120051 ; iPrescott et al.l 120071 ). and a majority of 
UV continuum and Balmer line emissions might be completely enshrouded by dust thus of 
course can't be corrected. There is no doubt that the above effect will be more severe for the 
central regions than for the entire galaxies, thus it could cause a higher aperture MIR-to-Ha 
luminosity ratio than the total value, in accordance with the discrepancies shown above. 

Dwarf galaxies were shown to be obviously deviated from the star-forming galaxies in 
the L(MIR)-L(Ha) correlations (see Figure Ufa) and (b)). This large deviation could be 
mainly caused by the aperture correction effect discussed abo ve, since many dwarf galax- 



ies h ave intensive star formation in their nuclear regions (e.g.. lGu et allbood : iLisker et al.l 
20061 ). This interpretation can also be proved by the aperture L(8/im[dust])-L(Ha) corre- 
lation shown in Figure ID^d), in which dwarf galaxies are only slightly deviated from the 
star-forming galaxies. Alternatively, this slight deviation could be explained by the fact that 
PAH and warm dust emissions are relatively weak in low metallicity environments (e.g. 



Madded 1200(1 lEngelbracht et al.ll2005l : IWu et all 12009 : iGalliano et al.ll2005l : IWu et al.l 120071 : 
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Engelbracht et al. 2008) 



5.2. Estimation of SFRs from MIR Luminosities 



Based on a sample of Spitzer FLS-field galaxies, IWu et al.l (120051 ) showed that both 
Spitzer 8/im(dust) and 24/nn luminosities could be used to calculate SFRs of normal star- 
forming galaxies, and they gave corresponding MIR-derived SFR formulae calibrated by 
either Ha or radio luminosities. However, the size of Wu05's sample is a little small (only 
have < 80 galaxies) and the MIR luminosities of their sample galaxies only cover a range 
of two orders of magnitude. In this work, we construct a large sample of 379 star-forming 
galaxies with the MIR luminosities range of three orders of magnitude. Besides correlating 
well with Ha luminosities, the MIR luminosities of our sample galaxies are also found to have 
tight correlations with TIR and FUV luminosities. These correlations allow us to estimate 
and compare SFRs from MIR luminosities calibrated by different SFR indicators. 

Although the correlations between L(MIR) with L(Ha) & L(TIR) are somewhat nonlin- 
ear, the nonlinearities of the above two correlations are generally less than 10%. Therefore, 
we can adopt linear correlations for calibrating SFRs (SFRs derived from non-linear corre- 
lations between L(MIR) with other luminosities also were illustrated in Figure [T3l . Based 
on the three linear correlations of L(MIR)-L(Ha), L(MIR)-L(TIR) and L(MIR)-L(FUV), we 
obtained that the SFR conversion factors of L(24//m) are 7.15x10 s , 7.79x10 s , and 7.04x10 s 
L , and those of L(8/im[dust]) are 1.58 xlO 8 , 1.70 xlO 8 , and 1.46 xlO 8 L , respectively. These 
factors are in good agreement with those of L(MIR)-L(Ha) and L(MIR)-L (radio) correlations 
by Wu05 with differences less than 10% (Figured^]). In fact, the differences in L(8/im[dust])- 
L(Ha) correlations between Wu05's sample and ours are mainly because of the galaxies with 
very low MIR luminosities. However, it must be noted that the MIR SFR formulae (espe- 
cially for 8/xm(dust)) derived above c an not be applied to dwarf galaxies (and other galaxies 
with low metallicity; IWu et al.l 120071 ) . since they deviate much from the linear correlations 
for normal star-forming galaxies. However, the deviation of dwarfs to normal star-forming 
galaxies in Figure [5] is apparently smaller than in Figure HI Therefore, it may be a simple 
and feasible method of using the combined observed Ha + MIR luminosities to compute 
SFRs. The black so lid line in Figure [TBI show the equation to compute SFR based on 24/zm 
luminosity given by lAlonso-Herrero et al.l (j2006al ) (their equation 3), which was derived from 
the correlation between 24/im and Ha luminosities for 30 L IRGs in the local universe, and 
this equation is similar to the one by ICalzetti et al.l (120071 ) except for different IMFs. We 
should keep in mind such large differences for infrared luminous galaxies. 



For TIR emission, the contribution from dust heating by the old stellar population 
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could not be ignored (~30% for normal galaxies, iBelll (120031 )); while for UV emission, except 
for the pollution from old star, the effect of intrinsic extinction is even stronger than that 
for Ha emission. Although the standard deviations of the fitting residuals from TIR (0.07 
for 24/xm, 0.06 for 8/xm(dust)) are smaller than from Ha (0.16, 0.18 respectively), the Ha 
emission trace the HII region surrounded massive star more directly. Therefore, we suggest 
that the MIR-derived SFR based on the MIR vs. Ha correlation is "better" than those 
derived from other correlations. 



5.3. Correlations among L(MIR), L(FIR) and L(TIR) 

The MIR luminosities L(MIR) are mainly from PAH and warm dust (VSGs) emissions, 
while FIR luminosities L(FIR) are dominated by cold dust emissions from large grains. From 
Figure El and comparing with the correlations between L(MIR) and L(Ha) shown in Figure H] 
and TableEl we find that the MIR lumin osities correla t e bett er with L(FIR) than L(Ha), and 



this result is consistent with that from iBoselli et all (120041 ) for late-type galaxies observed 
by ISO. 

However, correlations shown in Figures H] and [6] are on the basis of two different (MIR & 
FIR) samples. Therefore, we firstly checked the distributions of absolute B-band magnitude, 
redshift, u — r color, and EW(Ha) for both MIR and FIR sample galaxies (see Figure [2]), 
and found the two samples have similar distributions. To avoid possible selection effect, the 
L(MIR)-L(Ha) correlations for FIR sample galaxies are also shown in Figure [T41 The best fits 
for galaxies in the FIR sample are consistent with those for Wu05 and MIR sample galaxies. 
The standard deviations for the FIR sample are about 0.15, a little smaller than those for 
the MIR sample, since the FIR sample do not include galaxies with L(MIR) < 1O 8 L ; while 
their Spearman Rank-order correlation analysis coefficients are similar to those of the MIR 
sample. From the comparison described above, we confirm that the MIR luminosities are 
better correlated with FIR than Ha luminosities for star-forming galaxies. 

As described in §4.3 and shown in Figure among the four L(MIR)-L(FIR) correlations, 
L(24/xm)-L(70/xm) and L(8/xm[dust])-L(160/im) have much smaller scatters (the standard 
deviation of the fitting residuals are both 0.07) than the rest two correlations of L(24/xm)- 
L(160/im) and L(8/xm[dust])-L(70/im) (the standard deviation of the fitting residuals are 0.12 
and 0.10 respectively). Furthermore, the slopes of both L(24/xm)-L(70/xm) and L(8/xm[dust])- 
L(160/mi) correlations are close to unity, which indicates that such correlations are almost 
linear. The above results hint that the 24/im vs. 70/xm, and 8/im vs. 160/xm emissions may 
have similar physical origins: 

(a) The 24/zm emission is known to be dominated by warm dust emissions from VSGs, 
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which are mainly heated by young massive stars . Recent results by Spitzer ( jWu et al.ll2005 



Perez-Gonzalez et all 120061 ; ICalzetti et al.l 120071 ) have shown that the 24/im luminosity is 



one of the best SFR indicators. Previous works based on IRAS and ISO observations have 
proved that the FIR luminosity is also a good SFR tracer because of the emission peak 
around 60/im of dust heated by star formation, thus the 70/im emission must be closely 
related to star formation activities and should have tight and linear correlation with 24/im 
warm dust emission. 

(b) The 8/im ( dust) luminosity was also tho ught to be a SFR tracer by some authors (e.g., 
Wu et al.l 120051 ; lAlonso-Herrero et al.l l2006bl ). 8/im dust emissions include emissions from 
both PAHs and VSGs. In st ar-forming galaxies, the 7.7/im PAH emission dominates the 
8/im band (jSmith et al.1 120071 ) and its stren gth is likely to be affected by many factors, such 



2007 




Wu et al 


2007) 


(Peeters et al. 


2004 





as the radiation field and metallicity (e.g., iMaddenl |2000| ; iGalliano et al.ll2005t ICao &: Wu 



Moreover, PAH emissions trace B star better than young O stars 



tracer, and this could explain why the 70/im luminosities have better correlation with 24/im 
rather than 8/xm(dust) luminosities. 

(c) The 160/im emission is mainly from c old dust heated by old stars and is thought to 



locate in the diffuse regions of galaxies (e.g.. iGordon et al.1 120041 ); while the PAH molecules 



are generally located in the photo-dissociati on regions (PDRs), at the interface of HII regions 
(jHollenbach fc Tielens Ill997l ; lDraindl2003l ). The tight L(8/xm[dust])-L(160/xm) correlation 
indicates that the 8/im (dust) and 160/im emissions might be from the same regions of galax- 
ies, i.e., cold dust emission could exists in PDRs and it is also possible that there exists PAH 
emission i n the diffuse region s because PAHs can be pumped by cool stars with less UV 
radiation (ILi fc Draine II2002I ). Due to the low spatial resolution and shallowness of Spitzer 
160/im observations, firm conclusions on the physical connection between 8/im (dust) and 
160/im emissions and their distributions in galaxies must await high quality FIR data from 
future space telescopes (e.g., Herschel). 

Figure [6] shows that both MIR and FIR luminosities are well correlated with TIR lu- 
minosities for star-forming galaxies. Except for L(24/im), all of the 8/im (dust), 70/im and 
160/im luminosities are almost linearly correlated with the L(TIR). This result indicates that 
the monochromatic MIR and FIR lu minosities can be used for estimating the TIR luminosity. 



in ag reement with previous studies (jElbaz et al.ll2002l ; iTakeuchi et al.ll2005l ; iMarcillac et al. 
20061 ). 
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5.4. Properties of Composite galaxies and AGNs 



Investigating the properties of composite galaxies and AGNs will provide us some in- 
formation on the origin of MIR emission in AGN-hosting galaxies. Figures [8] and [10] show 
that composite galaxies and AGNs have distinct L(MIR)-L(Ha) and L(MIR)-L(FUV) cor- 
relations, i.e., they (especially for the ones with low Ha luminosity) have relatively higher 
MIR-to-Ha/MIR-to-FUV luminosity ratios than star-forming galaxies do. The high MIR- 
to-Ha/MIR-to-FUV luminosity ratios can not be caused by heavier dust obscuration in the 
host galaxies of composites and AGNs, since from Figure [3](c) it is clear that the host galaxies 
of composites and AGNs tend to be earlier Hubble types with redder u — r colors, and it 
is well-known that early- type galaxies should have less gas and dust (e.g., IWen et all 120071 ; 
Li et all 120071 ). Alternatively, for early- type galaxies in which the diffuse IR emissions can 
not be neglected, their MIR-to-Ha/MIR-to-FUV luminosity ratios would tend to be higher 
due to additional dust heating by abundant old stellar populat ions. Furthermore, the higher 
luminosity (Figure [3(a)) could re fer to higher me t allicityfe.g.. iLamareille et all 12004 ). while 
higher metallicity in early- types (IWu et all 120071 ; iLi et all 120071 ) could be another explana- 
tion to the higher MIR-to-Ha luminosity ratios of composite galaxies and AGNs except the 
re-emission of dust heated by AGN. 

From Figure [TT] it is clear that all the four L(MIR)-L(FIR) correlations for composite 
galaxies and AGNs are completely following those for star-forming galaxies. The L(MIR)- 
L(TIR) and L(FIR)-L(TIR) correlations for composite galaxies and AGNs are also shown to 
follow those for star-forming galaxies (see Figure [T2l) . Hence, the contribution from the active 
nuclear to both MIR and FIR emissions is more smaller than the contribution from the HII 
region, and it could be neglected for those weak AGNs. Nevertheless, the similarity of dust 
heating by AGN or HII region could be another possible explanation. These results indicate 
that we can estimate the total infrared luminosity accurately from monochromatic MIR and 
FIR luminosities, not only for star-forming galaxies (see also §5.3), but also for AGN-hosting 
galaxies. However, firm conclusion must await a quantitative analysis of a large, well-defined, 
and unbiased sample of FIR-selected galaxies from AKARI and Herschel observations. 



6. Summary 

We present and analyze the correlations between MIR, FIR, TIR, Ha, and FUV lumi- 
nosities for a large sample of galaxies selected from the Spitzer SWIRE fields. The main 
results described in this paper can be summarized as follows. 



1. The MIR (8/im[dust] and 24/im) luminosities of star-forming galaxies are found to be 
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well correlated with their Ha, FIR, TIR, and FUV luminosities. Based on these correlations, 
we derive corresponding formulae to calculate SFRs from L(MIR). All these formulae are 
well consistent with each other. 

2. The linear combination of observed Ha (Ha Q bs) and 24/im luminosities has tight 
and linear correlations with the extinction-corrected Ha luminosities L(Ha) for star- forming 
galaxies, so it can be used as a good SFR tracer. We determined the scaling parameter a 
in this combination, and also found the combined Ha Q b s + 8/im(dust) luminosities are also 
well correlated with L(Ha). 

3. Among the four L(MIR)-L(FIR) correlations for star-forming galaxies, the L(24/xm)- 
L(70/xm) and L(8/mi[dust])-L(160/zm) correlations are the tightest and almost linear. The 
former could be related to young massive star formations, while the latter might be relevant 
to diffuse dust emissions heated by old stellar populations. 

4. Composite galaxies and AGNs have higher MIR-to-Ha/MIR-to-FUV luminosity 
ratios than star-forming galaxies, nevertheless their correlations among MIR, FIR and TIR 
luminosities are completely following those for star-forming galaxies. 
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The BPT diag nostic diagram: [N il] /Ha vs [OIII]/H/5. The criteria from 



Kauffmann et al.l (120031 ) and lKewley et al.l (120011 ) are illustrated as dotted and dashed curves, 



respectively. The objects below the dotted curve are defined as star-forming galaxies, and 
the open circles represent dwarf galaxies with Mb > -18 mag. The boxes above the dashed 
curve denote AGNs, while the triangles between the above two curves are those classified as 
composite galaxies. 
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Fig. 2. — The distributions of (a) absolute B-band magnitude; (b) redshift; (c) u — r color; 
(d) EW(Hct) for galaxies in our MIR (black solid lines), FIR (red dashed lines), and UV 
(blue dotted lines) samples. The galaxy numbers of the above three samples are 605, 197 
and 421, respectively. 
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Fig. 3. — The distributions of (a) absolute B-band magnitude; (b) redshift; (c) u — r color; 
(d) EW(Ha) for 413 star-forming galaxies (black solid lines), 143 composite galaxies (red 
dashed lines), and 49 AGNs (blue dotted lines) in our MIR sample. 



-30- 




Fig. 4. — Correlations between MIR and Ha luminosities for star-forming galaxies. Panels 
(a) & (b): total Ha vs. total MIR luminosities; Panels (c) & (d): aperture Ha vs. aperture 
MIR luminosities. Dwarf galaxies are denoted as open circles. The best nonlinear and linear 
fits for normal star-forming galaxies (solid circles) are illustrated as solid and dotted lines. 
The dashed lines show in Panels (a) & (b) represent the nonlinear fits for FLS-field galaxies 
by Wu05, and the dashed lines shown in Panels (c) & (d) represent the nonlinear fits of total 
MIR-to-Ha luminositv correlations ( solid lines in Panels (a) & (b)). 
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Fig. 5. — Correlations between combined H« b s + 24/im & Ha Q bs + 8/xm luminosities and 
extinction-corrected Ha luminosities for star-forming galaxies. The symbols and line styles 
are the same as in Figure HI 




Fig. 6. — Correlations between MIR and FIR luminosities for star-forming galaxies. The 
symbols and line styles are the same as in Figure HI 




Fig. 7. — Correlations between MIR & FIR and TIR luminosities for star-forming galaxies. 
The symbols and line styles are the same as in Figure HI 
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Fig. 8. — Correlations between MIR and FUV luminosities for star-forming galaxies, com- 
posite galaxies and AGNs. The symbols are the same as in Figure [TJ The best nonlinear 
and linear fits for normal star-forming galaxies are illustrated as solid and dotted lines. 
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Fig. 9. — Correlations between combined FUV Q b s + 24/zm & FUV Q b s + 8/iin luminosities 
and extinction-corrected FUV luminosities for star-forming galaxies. The symbols and line 
styles are the same as in Figure [SJ 
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Fig. 10. — Correlations between MIR and Ha luminosities for composite galaxies (triangle) 
and AGNs (box). Panels (a) & (b): total MIR vs. total Ha luminosities; Panels (c) & (d): 
aperture MIR vs. aperture Ha luminosities. The nonlinear and linear fits for composite 
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Fig. 11. — Correlations between MIR and FIR luminosities for composite galaxies (triangle) 
and AGNs (box). The red solid and dotted lines represent the nonlinear and linear fits 
for composite galaxies, while the black dashed lines are the nonlinear fits for star-forming 
galaxies as shown in Figure [61 
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Fig. 12. — Correlations between MIR & FIR and TIR luminosities for composite galaxies 
(triangle) and AGNs (box). The red solid and dotted lines represent the best nonlinear 
and linear fits for composite galaxies, while the black dashed lines are the nonlinear fits for 
star-forming galaxies as shown in Figure [7J 
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Fig. 13. — SFRs derived from MIR luminosities based on the linear and non-linear corre- 
lations between MIR with Ha (red solid lines for linear, red solid curves for non-linear), 
FIR (green dotted lines for linear, green dotted curves for non-linear; notice that the lin- 
ear line and non-linear curve are coincident), and UV (blue short-dashed lines for linear, 
blue short-dashed curves for non-linear) luminosities. SFRs derived from MIR luminosities 
based on the relations from Wu05 between MIR with Ha (black long-dashed lines) and radio 
(black dashed-dotted lines) luminosities were also shown in this figure. The black solid curve 
in the left panel shows the e quation to compute SFR based on 24/xm luminosity given by 



Alonso-Herrero et al.l (l2006al ) (their equation 3). 
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Fig. 14. — Correlations between total MIR and total Ha luminosities for galaxies in the FIR 
sample. The symbols and line styles are the same as in Figure HI 
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Table 1. The expected and observed NIR — MIR colors of blue stars in aperture 

photometry 



Field 




Color 


Expected value 3. 


Observed value 


Corrected value 


EN1 


K 


- IRAC1 


0.04 


0.09 


0.05 


EN1 


K 


- IRAC2 


0.02 


-0.03 


-0.05 


EN1 


K 


- IRAC3 


0.00 


-0.03 


-0.03 


EN1 


K 


- IRACA 


0.00 


0.07 


0.07 


EN2 


K 


- IRAC1 


0.04 


0.06 


0.02 


EN2 


K 


- IRAC2 


0.02 


-0.03 


-0.05 


EN2 


K 


- IRAC3 


0.00 


-0.03 


-0.03 


EN2 


K 


- IRAC4 


0.00 


0.07 


0.07 


LH 


K 


- IRAC1 


0.04 


0.09 


0.05 


LH 


K 


- IRAC2 


0.02 


-0.02 


-0.04 


LH 


K 


- IRAC3 


0.00 


-0.02 


-0.02 


LH 


K 


- IRAC4: 


0.00 


0.08 


0.08 



from lEisenhardt et all (|2004l ): lLacv et all (|2005h 
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Table 2. The number of MIR, FIR and UV sample galaxies with different spectral types 



Sample 


Star — forming 


Composites 


AGN 


Total 


MIR 


413 


143 


49 


605 


FIR 


141 


49 


7 


197 


UV 


267 


65 


27 


359 



Table 3. Correlation coefficients for star-forming galaxies 



Sample 




y 


X 


a 


b 


s 


r 


c 


N 


(1) 




(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


FLS 


vL v 


24^to] (total) 


L{Ha](total) 


0.34±0.19 


1.12±0.07 




0.88 


1.29±0.19 


63 


FLS 


vL v 


8jj,m(dust)] (total) 


L[Ha](total) 


2.31±0.15 


0.92±0.05 




0.87 


1.68±0.19 


79 


MIR 


vL v 


24/xto] (total) 


L[Ha](total) 


-0.06±0.07 


1.18±0.02 


0.16 


0.92 


1.33±0.18 


379 


MIR 


vL v 


8jj,m(dust)] (total) 


L[Ha](total) 


1.14±0.07 


1.07±0.03 


0.18 


0.88 


1.67±0.19 


379 


MIR 


vL v 


24/um] (aper) 


L[Ha](aper) 


0.63±0.05 


1.10±0.02 


0.15 


0.93 


1.37±0.17 


379 


MIR 


vL v 


8fim(dust)] (aper) 


L[Ha](aper) 


1.74±0.05 


1.02±0.02 


0.14 


0.94 


1.88±0.16 


379 


MIR 


L[H 


a.obs /24^m] 


L[Ha](total) 


0.29±0.04 


0.96±0.01 


0.09 


0.95 


0.01±0.11 


379 


MIR 


L[H 


a>-obs/8fim} 


L[Ha](total) 


0.86±0.04 


0.89±0.01 


0.09 


0.94 


0.01±0.12 


379 


FIR 


vL v 


24^to] (total) 


L{Ha](total) 


0.28±0.10 


1.13±0.04 


0.15 


0.92 


1.37±0.17 


137 


FIR 


vL v 


8fj,m(dust)] (total) 


L[Ha}(total) 


1.87±0.11 


0.98±0.04 


0.15 


0.87 


1.69±0.17 


137 


FIR 


vL v 


24/im] (total) 


vL v \70[im\ 


-0.97±0.05 


1.04±0.01 


0.07 


0.98 


-0.56±0.08 


137 


FIR 


vL v 


8/j,m(dust)](total) 


vL„[70fim] 


0.87±0.07 


0.89±0.02 


0.10 


0.95 


-0.24±0.13 


137 


FIR 


vL v 


24/xto] (total) 


f-L„[160/zm] 


-2.75±0.10 


1.21±0.03 


0.12 


0.94 


-0.63±0.14 


137 


FIR 


vL v 


8jim(dust)] (total) 


vLv[l60fj,m] 


-0.69±0.06 


1.04±0.02 


0.07 


0.96 


-0.32±0.08 


137 


FIR 


vL v 


8[im(dust)](total) 


TIR[3-ll00^m] 


-0.44±0.05 


0.98±0.02 


0.06 


0.97 


-0.70±0.08 


137 


FIR 


vL v 


24/xto] (total) 


TIR[3 - 1100/im] 


-2.48±0.05 


1.14±0.02 


0.07 


0.98 


-1.02±0.09 


137 


FIR 


vL v 


70/im] 


TIR[3 - llOO/xm] 


-1.45±0.06 


1.09±0.02 


0.07 


0.98 


-0.45±0.09 


137 


FIR 


vL v 


160yum] 


TIR[3 - 1100/im] 


0.22±0.04 


0.94±0.01 


0.05 


0.99 


-0.38±0.06 


137 


UV 


vL v 


24/xto] (total) 


FUV 


-3.89±0.14 


1.29±0.04 


0.22 


0.86 


-0.91±0.24 


248 


uv 


vL v 


8[im(dust)] (total) 


FUV 


-2.58±0.14 


1.20±0.04 


0.22 


0.84 


-0.60±0.23 


248 



Note. — Col.(l): Sample used for correlation analysis. "FLS" represents the Spitzer First Look Survey sam- 
ple by Wu05; Col.(2)-(3): names of multi-A luminosities, total and aper represent the total and aperture lu- 
minosities, and all the luminosities are in units of L©; Col.(4)-(5): the coefficients a and b of the nonlinear fit: 
log 10 (y) =a+blog 10 (x); Col. (6): the standard deviation s of the fitting residuals; Col. (7): the coefficient c of the linear 



fit: log 10 (y) =c+log 10 (x); Col. (8): the coefficient r of the Spearman Rank-order correlation analysis; Col. (9): the 
number of sample galaxies (after excluding dwarf galaxies) used for the fitting procedures. 
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Table 4. Correlation coefficients for composite galaxies 



Sample 




y 


X 




a 






b 




s 


/• 




c 




N 


(1) 




(2) 


(3) 




(4) 






(5) 




(6) 


(7) 




(8) 




(9) 


MIR 


vL v 


'24/zm] (total) 


L[Ha](total) 


1 


95±0 


11 





96±0 


04 


0.19 


0.87 


1 


61±0 


22 


143 


MIR 


vL v 


8/im] (total) 


L[Ha] (total) 


3 


07±0 


12 





84±0 


04 


0.21 


0.79 


1 


86±0 


27 


143 


MIR 


vL v 


24fxm] (aper) 


L[Ha] (aper) 


1 


90±0 


10 





97±0 


04 


0.18 


0.89 


1 


64±0 


20 


143 


MIR 




8/im] (aper) 


L[Ha](aper) 


3 


11±0 


08 





85±0 


03 


0.15 


0.89 


2 


04±0 


21 


143 


FIR 


VLy 


'24/j.m] 


vL v \J0[im] 


-0 


75±0 


15 


1 


02±0 


05 


0.09 


0.90 


-0 


57±0 


10 


49 


FIR 


vL v 


8/tm] 


vL v [70fj,m] 


1 


70±0 


21 





80±0 


07 


0.14 


0.78 


-0 


36±0 


19 


49 


FIR 


vL u 


'24/zm] 


vL^[l60fim\ 


-2 


90±0 


31 


1 


22±0 


09 


0.16 


0.78 


-0 


60±0 


19 


49 


FIR 


VLy 


8/im] 


vL v [\Q§[im\ 


-0 


22±0 


1G 





98±0 


05 


0.09 


0.92 


-0 


38±0 


09 


49 


FIR 


vL v 


8/*m] 


TIR[3- 1100/im] 


-0 


02±0 


17 





93±0 


05 


0.10 


0.91 


-0 


79±0 


11 


49 


FIR 


vL v 


'24/tra] 


TIR[3- 1100/xm] 


-2 


69±0 


17 


1 


16±0 


05 


0.09 


0.90 


-1 


00±0 


11 


49 


FIR 


vL v 


70/im] 


TIR[i- 1100/xm] 


-2 


15±0 


15 


1 


16±0 


04 


0.08 


0.93 


-0 


43±0 


10 


49 


FIR 


vL v 


160/tm] 


TIR{3- 1100/xm] 





26±0 


12 





94±0 


04 


0.07 


0.95 


-0 


41±0 


09 


49 



Note. 



The definition of columns are the same as in Table [3] 



